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Spatial 't Hooft loops of strength k measure the qualitative change in the behaviour 
of electric colour flux in confined and deconfined phase of SU(N) gauge theory. 
They show an area law in the deconfined phase, known analytically to two loop 

C^") ' order with a "k-scaling" law k(N — k). In this paper we compute the 0(g 3 ) 

^ i correction to the tension. It is due to neutral gluon fields that get their mass 

OO i through interaction with the wall. The simple k-scaling is lost in cubic order. The 

Q\ • generic problem of non-convexity shows up in this order and the cure is provided. 

£\J ' The result for large N is explicitely given. We show that nonperturbative effects 

' appear at 0(g 5 ). 

y—i '. 

(N ■ 

| 1 Introduction 

The deconfining phase of QCD describes a plasma of in principle nearly free quarks and gluons. 
The colour electric flux of quarks and gluons is no longer confined in tubes, but is Debye 
screened. The pressure is then approximately equal to that of a Stefan-Boltzmann gas, with 
the number of internal degrees of freedom of gluons (N) and of quarks (Nf) appearing : 



, — i 

X 



PSB = (N*-l+ 7 -N f )^T* 

Of course there are interactions, becoming stronger as we go down in temperature. As the 
critical temperature T c is on the order of A-p-g ~ 200MeV the gauge coupling constant g 2 (T) is 
still 0(1) at T ~ 4T C which is where experiment may eventually take us. Unfortunately at these 
temperatures the perturbative pressure is wildly varying when adding higher order terms 1 2 . 
Since long one suspects the contributions from the long distance scales (Debye and magnetic 
screening lengths) to be the culprits for this failure of perturbation theory 3 . Lattice evidence 45 
suggests the series behaves well, once the order where magnetic scales appear has been taken 
into account. A prime example is the spatial Wilson loop: its dominant contribution is from 
magnetic scales, and it is accurately described down to T ~ 2T C 6 by letting the coupling run 
due to the modes of order T. 

Another example is the Debye mass. To lowest order it is given by excitations on the scale 
T. But already in next to leading order the magnetic scales contribute 8 . Once they are taken 
into account (and for any reasonable temperature they dominate 9 ) the remnant of the series 
is small according to numerical simulations 4 5 . 

The magnetic scales contribute a partial pressure of order (g 2 T) 3 , so come in only at 0(g e ) 
in the pressure 10 . They give a non-perturbative contribution to this coefficient, which is under 
study 73 . 



1 



So the data seem to tell us that the magnetic sector, in establishing a starting point for 
the series, is just as important as the electric sector with its gluonic quasi-particles. 

A rationale may be provided by simply postulating magnetic quasi-particles 11 with a 
density ~ (g 2 T) 3 . A simple consequence is k-scaling for spatial Wilson loops at very high T. 
This was confirmed by lattice data 12 to one percent accuracy. 

In this context it is of interest to compute analytically the low orders of the spatial 't Hooft 
loop, to compare the result to lattice simulations, so as to see from what order on the series 
behaves well. 

This loop measures the colour electric flux going through it. One can write it as an electric 
dipole sheet in say an x-y plane at some fixed z and t coordinate. On the lattice the dipolc 
sheet appears as an electric twist on that same x-y cross section of the periodic box. 

The pressure is not changing, but a boundary effect occurs, proportional to the area of the 
twist. The effect of the electric dipole sheet will extend over a distance typically of the electric 
screening. So it is reasonable to expect that the magnetic scales, much larger then the Debye 
screening, will appear only through their effect on the Debye screening (which at T <~ 2T C is 
dominating the lowest order Debye screening 9 ). We find that this effect comes in at 0(g 5 ), see 
section(4). 

Hence we expect to produce a purely perturbative series up and including 0(g 4 ). Numerical 
simulation is available 18 from very near T c , where perturbative methods are certain to fail, to 
about 2.5T C , where the data are in reasonable agreement with two loop results 14 . In order to 
do a more meaningful comparison lattice simulations of the loop from 2.5T C on should become 
available. 

Also in this paper we extend the computation of the 't Hooft loop from the known two loop 
order g 2 to include all g 3 effects. Apart from the motivation above, this computation serves a 
two-fold purpose. First it will verify or falsify the simple k-scaling in the strength of the loop, 
found to one and two loop order. 

Second there is a general motivation. The computation of the tension involves a potential 
between two degenerate minima, hence a potential with a non-convex part. If the strength of 
the loop is k the profile induces masses of the order O(l) in the coupling in 2k(N — k) gluon 
fields. The others stay massless, since they do not interact with the wall to one and two loop 
order, so do not contribute to the tension. This is the cause of "k-scaling" . But in three loop 
order they can develop a selfenergy, through which they start to interact with the wall. The 
order g 3 involves the the resummation of the Coulomb propagator by this selfenergy. The latter 
is related to the second derivative of the potential, hence with a negative part. The strength of 
the loop enters in this selfenergy and it turns out that to leading order in the number of colours 
there is a window in the strength where the selfenergy stays non-negative and where we can 
solve for the tension. That is what we will do in this paper. 

On the other hand the generic solution to the problem is of course suggested by the physics 
of the problem: the selfenergy is a consequence of the collective interaction of the wall with 
the massless gluonficld. Hence we have to substitute the profile into the selfenergy and solve 
for the eigenvalue spectrum of the Schroedinger equation. It turns out this spectrum can be 
found with the methods of reference 23 . Hence it is interesting to see how the two approaches 
compare in their common domain of validity. This will be done in a sequel to this paper 24 . 

The lay-out of the paper is as follows. Section(2) sets the problem in the familiar context 
of a spontaneously broken centergroup symmetry. Then in the next section the definition of 
the 't Hooft loop is given, and an intuitive argument for its behaviour in hadron and plasma 
phase. 

Then we start the main thrust of the paper in section(4) with an outline of our semiclassical 
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approach and the necessary tools. One, and two loop results are reviewed in the next sections. 
The infrared divergencies in the three loop case are discussed and how to get rid of them 
respecting centergroup invariance. Then in section(6) we compute the 0(g 3 ) terms for large N 
and a window of 't Hooft loop strengths. 

In section(8) the results are discussed. Appendices contain technical details. 



2 Z(N) symmetry in hot gluodynamics 

In absence of matter fields transforming non-trivially under the centergroup Z(N) of SU(N) 
there are two distinct canonical symmetries: magnetic and electric Z(N) symmetry 13 17 . They 
act through charge operators on the quantum states of the theory. 

There is another Z(N) symmetry 13 induced by a compact periodic dimension. Only the 
timelike case is here of interest. It is a symmetry of the pathintegral giving the pressure of 
the theory. This path integral is periodic in Euclidean time. It acts on the fields of the path 
integral as a gauge transformation Q[ with a discontinuity expiZ^ in the centergroup. Since 
the fields do by assumption not feel the action of the centergroup their periodicity is respected 
by this gauge transformation. The action is invariant and so is the integral. To make this into 
something more useful we have to have an order parameter, that does transform non-trivially. 

This is the timclike Polyakov loop: 

P(A ) = 1-TrTcxpi [ drA ( T ,x) (1) 



It transforms in an obvious way: 



P(i4°')=exp«^P(4>) (2) 



So the constrained path integral we get by integrating over all configurations with a fixed value 
of the Polyakov loop will then have the same value in all Z(N) transformed points, because only 
the loop is transforming. Define the effective action, with matter fields represented by dots: 

ex p-l[/(P) = J DA...S(P-P(Ao))e3qp~S(A,...) (3) 

then 

U(P) = U(expil^P) (4) 

The effective action is Z(N) invariant. P is the VEV of the Polyakov loop. If the VEV is 
non zero, like in the deconfined high T phase, we have a Z(N) periodic potential for the VEV. 

This potential is quite useful for the physics of the plasma. First of all it can be seen on 
general grounds, that the potential must have a minimum for the space averaged loop P = l 21 . 
The potential is known to two loop, 0(g 2 ), and we will compute it below to 0(g 3 ). It is gauge 
invariant by construction and may serve as a definition of the free energy at P = 1. This 
definition has the advantage that the Feynman rules have a built-in gauge invariant infrared 
cut-off. It is given by by the scale qT, which determines the deviation of P from its bulkvalue 
1: P~ exp-qT/T. 
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Figure 1: The 't Hooft loop in an x-y cross-section of the box. Arrows indicate the electric fieldstrcngth projected 
on the hypcrcharge Yfc. Eventually we will consider the loop stretched to the full cross-section. 



3 The 't Hooft loop and colour electric flux 

The behaviour of colour electric flux in the confined and deconfined phase is expected to be 
qualitatively very different. In the confined phase we expect rings of electric fluxtubes to 
be the glueballs, whereas in the deconfined phase the gluons are the electric quasi-particles: 
approximately a free Bose gas. 

The 't Hooft loop Vk(L), the closed loop C say being square, size LxL in the x-y plane, 
is usually defined 13 as a closed magnetic flux loop L of strength k^, with N the number 
of colours. In operator language it is defined as a gauge transformation with a discontinuity 
expikjf when crossing the minimal surface. As we will justify below, one can write this gauge 
transform as a dipole sheet 1722 , see fig.(l): 



with N-k entries k and k entries k-N, to have a traceless matrix. The charges Yfc are gener- 
alizations of the familiar hypcrcharge, with k=l. The charge Yfc of a gluon is or ±N. The 
multiplicity of the value N is k(N — k). The same is true for the value — N. So, e.g. for N = 3 
and k = 1 one finds the four kaons with hypercharge ±3. 

Exponentiation of Yfc gives expi^Yfc = expik^ = Zk, the centergroup element. 

E z is the z component of the canonical electric field strength operator E = \ a E a . 

The 't Hooft loop Vfc(L) can be written as expz^^fc with 




(5) 



where the NxN diagonal traceless matrix Yfc is defined as 



Yfc = diag(k, k...k, k-N,k- N , fc - JV) 



(6) 




S(L) 




(7) 



the electric flux operator. 

If a spatial Wilsonloop in the fundamental representation - a fundamental electric flux loop 
- W(L') = expi J L , dl.A is piercing through the area subtended by the 't Hooft loop, the loops 



"The A matrices being normalized to Tr\ a \b 
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do not commute anymore with each other. The effect'' is a phase expik^- and can be written 
as the 't Hooft commutation relation: 

Vki^WiL^VkiL)- 1 = expz— Y k W(L') (8) 

The appearance of the phase in the commutation relation shows that the 't Hooft loop is 
a gauge transform with a discontinuity expi^, detected by the Wilsonloop. 

The choice of the charge Yk is determined up to a regular gauge transformation. That 
will not change the commutation relation above. In particular a permutation of the diagonal 
elements will leave Zk invariant. The reader might object that to a given centergroup element 
there corresponds a lattice of points in Cartan space. So which point to choose? The answer is: 
the assignment of a specific value of Zk to the plasma groundstate is arbitrary. But the jump 
at the dipole layer must be such that one arrives at a nearest neighbour of the chosen lattice 
point c . 



3. 1 Thermal expectation value of the 't Hooft loop in the quasiparticle picture 

Let us imagine the hadronphase at some temperature T - below T c - as a gas of microscopic 
closed fundamental flux loops, i.e. Wilsonloops on the scale of glueballs. Then, from the 
commutation relation above, we see that only glueballs sitting on the perimeter L of the 't 
Hooft loop will have an effect on the 't Hooft loop. 

But, on the contrary, a gas of free gluons will change the behaviour of the 't Hooft loop 
into an area law behaviour. With the tension p k (T) and the area A(L): 

<V k (L) >= cxp - Pk (T)A(L) (9) 

This area law follows from the quasi-particle picture in the following simple way. 

Just one single gluon will change the value of the loop by a minus sign. To see this, note 
that the gluon has a charge or N. This is true for any value of k in view of the definition 
of k charge, eq.(6). What changes is the multiplicity of the charge ±N. We know this charge 
is screened by the plasma, with a screening length Ie- So if the particle is within distance Ie 
from the surface spanned by the loop one half of the flux will go through the loop, the other 
half docs not. As a result the loop captures a flux, eq.(7): 

= y (io) 

and the loop Vk(L) = expi^$fc acquires the value —1. 

Let P{1) be the probability that I gluons are in the slab of thickness Ie- We suppose this 
probability is centered around the average (l), with width proportional to (I), as behooves a 
thermodynamic distribution for a gas of bosons with a mass. The outcome for the average is 
then: 

^(-yP(Z) = exp-c</> (11) 

with the proportionality factor c related to the shape of the distribution d . This result is true 
for one gluon species with charge ±N. We know from eq.(6) there are 2k(N—k) of such species. 

^through the canonical commutator [E?(x), A* (y)] = i8 a ' b 8^ n 8(x — y). 
c An intriguing explanation has been given by Kovncr 20 . 

d The distribution for massless bosons does not have this shape and does not give an arealaw. 
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Since they are acting independently the tension becomes: 



Pk = c2k(N - k)l E n{T) 



(12) 



The dependence on the strength k is the k-scaling law for 't Hooft loop tension. The 
screening length is proportional to -, 1 . 



Summing up: a gas of free gluons predicts the ratio Pk/pi = k(N — k)/(N — 1). The factor 
k(N — k) is just the number of gluons that have Yfc charge N. 

4 Semi-classical determination of the 't Hooft loop expectation value 

The arguments in the previous section were semi-classical, in that only the density n(T) of the 
free gluons is quantum mechanical. So we can expect that a semiclassical approach would be 
in place with a systematic expansion in the coupling. 

Apart from its obvious connection with colour electric flux the 't Hooft loop is intimately 
related to the Z(N) symmetry of section(2). Take its order parameter P(Ao), an d move it 
through the dipole layer. Then it will get multiplied by the discontinuity expik^, being a 
(heavy) fundamental test charge. Immersing the loop in the plasma induces a disturbance. 
The disturbance is described by a profile C. Since the loop is gauge invariant the response of 
the plasma is too. This profile is the phase of the Polyakov loop as a function of its distance 
to the minimal surface. 

So the approach will be to compute the free energy excess AF due to the presence of the 
Polyakov loop profile. Let the box be of size L\ r xL z , with L z » L tr and both macroscopic. 
Extend the loop in fig.(l) to the full x-y cross-section of size L\ r and located at say z = 0. 

Then we have exp— AF(C) — J DC exp — -^fU{C). U{C) is the constrained effective 
potential we mentioned in section(2). 

We will minimize this excess free energy by varying the profile C, find the profile C m in and 
the minimum free energy U(C m in), which is the tension of the 't Hooft loop. 

This is a simple procedure. What is technically involved is the transformation of the gauge 
potential to the gauge invariant Polyakov loop. This is done in the next subsection. 

4-1 The constrained pathintegral 

So the thermal average of our loop is given by the Gibbs sum 



The average can be transformed into Euclidean path integral form. When doing so, one has 
to realize that the operator expi^ J s (L)dxdyTrE z Yk is only gauge invariant on the physical 
subspace. This follows from it being a gauge transform with a discontinuity in the center of 
SU(N). The discontinuity commutes with all of the gauge group. If we chop up the operator 
in little Euclidean time bits St, then every bit apart will not be a centergroup element, hence 
will not be gauge invariant!. This in contrast to the Boltzmann factor where this operation will 
not affect the gauge invariance. So the operator stays inserted at a fixed time slice, say r = 0. 
Only there Aq ^ 0, for all other times A$ = 0. We get: 




< Vfc(L) >= TrV k (L)exp~H/T 



(13) 



< 




DAoDAexp-—^ 
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The terms in the exponential constitute the action S(A, Sk) with the original surface dis- 
tribution of dipoles as a source sf = 2irYkdi6(z) localized at z = 0. This source induces a 
profile in the thermal Polyakov loop, or more precisely in the average over transverse directions 
P(A ) = ^JdxdyP(A ). 

The source term tells the Polyakov loop to jump over P(^Yk) — expik^- at the minimal 
surface and to have value 1 at long (longer than the Debye length Ie distance from the surface. 

So the thermal average becomes in terms of the profile P(z): 

<V k (L) >= J DP(z)exp-^rU(P) 
The effective action U(P) is given by the constrained path integral we introduced in section(2): 

cxp-^|[/(P) = J DA DA6(P-P(Ao))e X p^S(A,s k ) (16) 

and we have used an abbreviated notation for the constraint: 

S(P - P(A )) = n Zil 6(PV\z) - P l (A (z))) (17) 

where / runs from 1 to AT — 1. In the following we parametrize the fixed loop by a diagonal 
traceless NxN matrix C{z) = diag(Ci, C*2, Cat) : 

= l-TrexpilC(z) (18) 

So the effective potential is defined on the Cartan space in which the matrices C live. The 
matrix describes the profile of the loop. Once again, the way the strength k comes into the 
effective action is through the boundary condition at the dipole layer. 



(15) 



4.2 Thermodynamic limit and minimizing path 

In the thermodynamic limit L z and the transverse size L tr become very large with respect to 
any microscopic scale. 

Then the integral over profiles, eq (15), gets its main contribution from the minimizing 
configuration C m i n . Writing the fluctuations as C(z) — C m i n (z) + 7(2) we obtain: 

< Vk(L) >= exp-^U(C min ) J ^ 7 ( x )exp™^| 7 W^"(C mm ).7(^) + 0((7) 3 ) (19) 
So the tension will become: 

Pk = U(C min ) + -LTrlogU"(C mm )) (20) 

The stability matrix U" has non-negative eigenvalues for all fluctuations in the profile, because 
our kink C m i n is stable e . 

The strength of the loop determines the jump at the dipole layer. It also fixes the path in 
Cartan space (the space of diagonal C's) along which the minimal profile is realized. The path 

e No zeromode appears because the profile location is fixed by the fixed loop. 
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turns out to be the simplest possible n : if parametrized by q, < q < 1, it is given by the one 
dimensional set of Cartan matrices 

Y k (q)=qY k (21) 

with Yk the charge characterizing the strength of the dipole layer, eq.(5). In exponentiated form 
it goes from 1 to expi/c^, as q goes from to 1. Z(N) invariance of the profile functional U(C) 
garantuees we we can take a smooth path from 1 to expik^-, instead of a path that makes a 
jump expik^- at the loop and returns to 1. A proof of the rectilinear path being the minimal 
one is still lacking. Only for SU(3) and SU(4) it is known to be the case by inspection 14 11 and 
a proof at large N is given in ref. 14 . 



4-3 The gradient expansion of the effective action 

In principle the tension is given by eq.(20). We have now to work out an approximation scheme 
in weak coupling. 

The effective action appears in the constrained path integral. This path integral contains 
the profile P(z) as parameter. As long as we limit ourselves to profiles with gradients much 
smaller than the profile, a gradient expansion in d z P(z) should make sense. 

To see this go from the A a = gauge in eq.(14) to background gauge. So the potentials 
and therefore the Polyakov loop are parametrized through a colour diagonal background field: 

A fl = B(z)S IMfi +gQ ft (r,x) (22) 

and choose background gaugefixing Tr{D(B)^Q^) 2 . A systematic loop expansion of the con- 
strained path integral is possible 19 . 

To lowest order the background field B equals the phase C of the Polyakov loop profile 
through the constraint in the pathintcgral, eq.(16). The phase C appears in all derivatives 
Dq(C). Since the phase is diagonal in colour the colour basis for the quantum fluctuations 
in eq.(22) that diagonalizcs the covariant derivative is the Cartan basis: A y are the N(N-l) 
matrices with only one entry (i ^ j) non-zero and equal to Then the (N-l) diagonal 

matrices X d with the first d-1 entries equal to 1 , the d'th entry equal to 1-d, and all others 
equal to 0. The normalization stays \ as for the Gell-Mann matrices. The fluctuating fields Q 
can be decomposed on this basis: 

ij d 

The covariant derivative will act as Dq{C) = do + i(Ci — Cj) on the components and like 
Dq(C) = do on the diagonal components. So the Feynman rules are the usual background field 
rules, with the Matsubara frequency po shifted by the phase of the Polyakov loop. The phase 
C{z) commutes with d z in the gradient expansion. For the Coulomb part of the propagator we 
use the resummed propagator: (Aqq) -1 = p 2 +m 2 D + T 2 (d — Cj) 2 for the Q y propagators and 
(Aq ) _1 = p 2 + m 2 D for the Q d propagators. 

We are interested in the effective potential along the minimizing path Yk{q). Then a 
simplification occurs in the Feynman rules. Along that path we have Dq(C) = do ± i2wq for 
those gluon fields Qij with charge ±N. If the charge is zero the covariant derivative reduces to 
do- So in fig. (2), (3), (4) and (5) the lines and vertices obey the above rules. 
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c) 



d) 



Figure 2: (a) is the one loop contribution, and stands for gluon and ghost loop, (b) are the two loop contributions, 
(c) is the rcnormalization of the Polyakov loop (double circle) inserted (dotted line) into the one loop, see text 
below eq.(34). (d) is the kinetic term with the loop being a gluon or ghost loop 



4-4 One loop tension 

The only classical term that survives in the effective action is the electric field strength: 



No potential term survives in U c i(C). And this means there is no area law in this classical 
approximation: in our box of size L z L\ r the profile minimizing U c i is linear in z. Hence the 
profile yields a flat action density, so no wall results. 

So we have to expand to one loop order from the graph in fig. (2a). 

The gradient expansion gives a term proportional to Tr(d z C) 2 , and a potential term with- 
out gradients: 



To find the minimum configuration for this form of action is well known 14 . The result is: 



U(C) = U cl (C) = ^Tr(d z C) 2 



(24) 




(25) 



-d z C(K(C))^-(V(C))i=0 



(26) 



as equation of motion for C = C\ 
The tension reads: 




(27) 



We now write 



K = 



V = 



l + g 2 K ± + 
V 1 +g 2 V 2 + 



(28) 
(29) 
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In the one loop approximation the potential term equals: 



Vi(C) = ^J2(^ 2 ( 1 -\^\) 2 ( 3 °) 

ij 

where dj — d — Cj . 

Now we focus on the minimizing path Yk(q). Then, from the explicit form of the charge 
Yfe, eq.(6), we see that all non-zero dj — ±2irq. As the effective potential V\(C) is even in dj 
all terms contribute equally. And there are 2k(N — k) pairs ij with non-zero charge. Similar 
counting applies to the classical kinetic term, and after these manipulations we have for the 
sum of the classical kinetic term and the one loop potential V\(C) on the path Yk(q): 

rp2 f 1 97T 2 T 4 

-Tr(^C) 2 + V^C) = 2k(N k) J dz( — (2nT) 2 Tr(d z q) 2 + -^ 2 (1 - kl) 2 ) (31) 

So we have found once more the k-scaling of the quasi-particle picture in section(2)! The 
quasi-particles with charge ±N correspond precisely to the fields with massive propagators, the 
"neutral" particles to the massless propagators. Only massive propagators contribute to this 
order. 

Searching for a minimizing profile q m in(z) forces us to to choose a configuration q m i n (z) 
with d z q m i n (z) — 0(g). This in order for the electric term to be of the same order as the 
one-loop term. So it vindicates the use of the gradient expansion. The dependence on z of the 
profile is through the combination gz. Of course this was to be expected, because the electric 
dipole layer disturbs the plasma only over the Debye screening distance. 

Still we did not mention the boundary condition at the surface spanned by the loop. As 
in a surface layer of dipolcs in electrostatic Maxwell theory we expect the electric field to be 
continuous when going through the layer. This then fixes the minimal profile. 

The equations of motion for the profile become 

d z q - m D q(l - q) = (32) 

for q = q min with m 2 D = ^-T 2 the Debye mass. 

Note that the equations of motion and hence the profile are independent of the strength k. 
The tension of the wall in this leading order is given by 14 11 

= k(N - k) 4?f2 J 2 (33) 
k 3(3 5 2 A)5 

The tension has the same parametric dependence on N, k and the 't Hooft coupling g 2 N 
as we found in the quasiparticle picture. 



4-5 Next to leading order contribution to the tension 

The correction to the leading result comes from combining the one loop kinetic term with the 
two loop potential term evaluated in the background q(z). This will give an 0(g 2 ) correction 
to the result above and has been discussed elsewhere 14 19 11 . Here we only look in some detail 
to the infrared aspects. 

To begin, the two loop correction V^c?) to the potential is infrared finite and equals (see 
fig.(2b)) 

9 2V2 W = - 5 j£v k ( N fc )« 2 (! - M) 2 (34) 
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Important is to realize that the Fcynman rules from the constrained effective action eq.(16) are 
not only involving the traditional vertices, but extra vertices due to the contraint P — P(Aq). 
They come on one hand from expanding -P(A)) in the quantum field, and on the other hand 
from insertions into traditional vertices. They produce the effects of the renormalization of the 
Polyakov loop. In the two loop approximation they give rise to only one diagram, representing 
the effect of the one loop renormalization of the Polyakov loop inserted in the one loop free 
energy, fig. (2c). 

The result from fig. (2c) was discussed in ref. 14 and in a systematic way in ref. 19 . As far 
as the tension concerns there is no qualitative difference with the leading term, and the result 
turns out to be again k-scaling 14 11 (see appendix A): 

Pk {T) = pi 1} (T)(l - ^iV (15.2785... - f( lE + 1)) + 0(g 3 ) (35) 

The coupling is now running, due to the one loop kinetic term. It is defined in the MS scheme 3 . 



4-6 Profile to two loop order and Debye mass 

The behaviour of the corrected profile far from the loop is less straightforward. This is obvious 
since the Dcbye mass controls the behaviour there, and the next to leading correction to the 
Debye mass gets contributions from the electrostatic sector, and from an infinite number of 
diagrams in the magnetostatic sector 10 . 

So let us look at the contribution from the static sector to the kinetic term K i (q) due to the 
diagram in fig. (2d). From the Feynman rules in Feynman background gauge we find the only 
contribution from a Coulomb and a spatial propagator, carrying mass m 2 D + (2nTq) 2 7 (2irTq) 2 
respectively: 

<?Wi(<z) = Ag 2 NfT [ U (36) 

J (2^) 3 (l 2 + m 2 D + {2^Tqf){{l+p)^ + {2wTq) 2 ) 2 

The external momentum p = \p\ is O(gT) as we argued above for the gradient of the profile. 
In this section the mass parameter q — 0{g 2 ) will be parametrically smaller than the Debye 
mass, hence we will neglect it in the Coulomb propagator. After integration one finds: 

f Kl (q) = -4f^(llog - iP + m ° + 2 ^J ) (37) 
y W 1 8tt y ip y ip + m D + 2irqT J V ' 

where we put the cuts along the imaginary axis in the complex p-plane, starting from ±.i(mu + 
2tt qT). 

So the pole ip — mo is well separated from the cut in eq.(37), so we substitute mj for 
ip. The gradient expansion for the static modes has changed the 0(g 2 ) graph into a term 
linear in the coupling g, multiplied by the logarithm. Put into the equation of motion for the 
profile, eq.(26), it gives a correction to the Debye mass appearing in the lowest order equation 
of motion, eq.(32). When q = 0(g 2 ) we get Rebhan's correction 8 to the Debye mass: 

[W q 2 N 1 
m D = J-gT+ y —log{-) (38) 
V o 47r g 

There is also a static contribution in fig. (2d) from the spatial propagators only. It reads 

g 2 fK 1 (q) = -\g 2 NfT [ i (39) 

3 J (2n) 3 (P + (2irTq) 2 )((l+p) 2 + (2nTq) 2 ) 2 
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Doing the integral will give us an amplitude with cuts starting at ±iAirqT and our procedure 
will fail for q on the order of the magnetic mass because the pole is now on the cut. This failure 
indicates that other means than perturbation theory are necessary 10 . 

It should be emphasized that only when q = 0(g 2 ) we get this non-perturbative be- 
haviour 10 . For q parametrically larger the result is still perturbative. So in the q-integration 
leading to the tension, eq.(27): 



the non-perturbative contribution will contribute to 0(g 5 ). 

In conclusion, the next to leading order has electrostatic and magnetostatic divergencies 
that show up only in the profile at long distance, m£,\z\ — log(p-), and in the tension at 0(g 5 ). 
They typically show up for profiles on the order of q — 0(g 2 ). At values of q — 0(1) the 
theory is protected from the infrared, to this order. This concludes the discussion of the next 
to leading order diagrams. 

5 Linear divergencies in three loops and Z(N) invariant selfenergy 

In this section we analyze the two loop kinetic term and the three loop potential term. We 
look for infrared divergencies that might contribute to cubic order to the tension. It will turn 
out that none of those are present in the kinetic term, but indeed in the potential term. 

5.1 Three loop potential and self energy matrix 

From the preceding section one would conclude that only the fringes of the wall, near the 
unperturbed plasma, have infrared divergencies. The three loop potential, surprisingly, shows 
them for all values of the profile-^. Let us look in more detail to this. 

The three loop diagrams are divided in two sets, Df and D p . The diagrams in Df are all 
diagrams with the topology of free energy diagrams, but put in the background q of the tunnel- 
ing path. Their calculation involves only vertices from the Yang-Mills action, background gauge 
condition and ghost action. D p involves the vertices due to the constraint and renormalizes the 
Polyakov loop. 

(2) 

The set Df is divided in turn into two-particle irreducible diagrams Dy and one-particle 

irreducible diagrams D^p . The two particle irreducible diagrams have been analyzed for q = 
and are individually infrared convergent. This means they are 0(q) for small q, like their two 
loop counterparts. There is a single one-particle irreducible diagram as shown in fig. (3). 

This diagram is in the case of vanishing background linearly divergent for the 00 ( "Coulom- 
bic") components of the propagator. Similarly, the diagram in fig. (4) showing the self energy 
insertion on the one loop renormalization of the Polyakov loop has a linear divergence, which 
in the case of vanishing background has been analyzed by Gava et al. 15 . 

For notational convenience we write IL^.,(q) for the Coulombic selfenergy matrix at zero 
momentum and frequency, dropping the Lorentz indices. The index k means we are evaluating 
the selfenergy along the tunneling path Yk(q). If the two propagators have no background 
induced mass the diagram is linearly divergent in the infrared and we have to resum the prop- 
agator by including the selfenergy 11^, (q) at zero Matsubara frequency and momentum. At 

^In a previous paper 11 we overlooked this possibility 
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Figure 3: The only three loop diagram of free energy topology with an infrared divergence. The shaded blob 
is the one loop selfenergy of fig. (5). The colour index of the two propagators need not be the same due to 
background dependence inside the blobs. 



o 



Figure 4: Only three loop insertion diagram with infrared divergence 

q = we have the Debye mass 11^,(0) = m 2 D 8 CyC i. So the resummation we used up to now 
subtracts correctly the selfenergy at q = 0. But to have a finite result for q ^ we have to 
subtract the selfenergy for all q, and hence resum it for all q. Only then the linear divergence 
is cancelled. So the inverse Coulomb propagator in the background q becomes: 

iXc + ngfa) (4i) 

and gives an order g 3 contribution to the effective potential. In particular for q = 0, this 
contribution is proportional to m D . 

For the diagrams in fig. (3) and fig. (4) to be finite we have to check that the selfenergy at 
small momenta behaves like p 2 , and this is indeed the case. 

For the infrared finite result to be Z(N) invariant we should keep contributions from all 
Matsubara frequencies in the selfenergy. After doing so the selfenergy reads: 

n$(«) - 4g 2 f c ' a ' b B 2 (C a )r'' a ' b (42) 

where .62(g) 1S related to the Bernoulli polynomial: 

B2{q) = ^(\-\q\+q 2 )^B 2 (q) (43) 

This Bernoulli polynomial is periodic mod 1, because we summed over all Matsubara frequencies 
in the selfenergy. The formula is only true for external legs c and c' being neutral, i.e. not 
feeling the background field. In that case C a = —Cb- It is proven in appendix B. Note that the 
dependence on the channel k comes in only through the argument C a - 

The argument of B 2 (C a ) is if the index a is that of a diagonal basis element, or d — Cj 
if a = (ij). The group structure constants are defined as: 

^a.b.c = 2Tr[\ a ,X b ]\ c 
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Figure 5: The (00) component of the gluon selfenergy of with fixed colours a and b on the internal lines as in 
eq.(42). See Appendix C for evaluation. 



Eq.(42) is diagonalized for any tunneling path in appendix C. 

Here we give the results for the eigenvalues 11^ (g). Obviously all eigenvalues have the 
Debye mass at q = and q = 1 in common. The selfenergy matrix contains in every channel k 
the N — 1 by N — 1 matrix of diagonal gluons. Furthermore the indices c and c' can take any 
value ij for which the background = in the incoming propagator. Thus the dimensionality 
of the matrix depends on k through the number of zero eigenvalues of the adjoint representation 
of Y k : N 2 — I — 2k(N-k). 

In terms of B 2 (q) = B 2 (q) — g and dropping the overall scale factor m 2 D one finds 
(k — l)(fc + 1)) eigenvalues 

(N-k-l)(N-k+l) eigenvalues 

and one eigenvalue 



N — k ~ 

(l + Q^ r )B 2 (q) (44) 



l + 6^B 2 (q) (45) 



l + 6B 2 (q) (46) 
5.2 Infrared behaviour of the Polyakov loop 

The graph in fig. (4) shows the renormalization of the Polyakov loop to 0(g 3 ). The calculation 
parallels that in ref. 14 , 19 , but now with the resummed gluon propagator from eq. (41), and 
the diagonalized version of eq.(42) as described above. The effect on the eigenvalues of the 
Polyakov loop matrix is given by: 



2tt (4tt) 2V s/ v 2tt ; (8tt) T 

The second term is due to heavy modes of order T and contains the gauge dependence 
needed to cancel the gauge dependence in the graphs of fig. (2). The cubic term 9 contains as 
first term the one discussed by Jengo and Gava 15 without background field. The second term 
P3 is induced by the background. The precise dependence on the background we do know but 
it is irrelevant for the purpose of this paper as the following argument shows. 

The cubic term is inserted into the lowest order result eq. (30) and gives a cubic contribution 
to the tension: SCijB^Cij). This cubic contribution, unlike the quadratic one, is zero! 
The reason is that the background contribution is only present in those combinations (ij) with 
Cij = 0, and -83(0) = 0. For the background independent correction the sum B 3 (dj) in 
the insertion cancels out on the path qY^. This is because k(N — k) terms multiply -83(g) and 
an equal number B 3 (—q). Since B 3 is an odd function 14 the result cancels out. 

'Remember m D /T = (S^L) 1 / 2 . 
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5.3 Infrared behaviour of the kinetic term to two loop order 

We still have to understand wether the two loop kinetic term can contribute to cubic order. This 
term, nominally of 0(g 4 ), will contain contributions of order g 3 . The question is wether they 
are present for all values of q including q = 0(1). So we will insist on at least one propagator 
with a background induced mass of order T, i.e. q — 0(1). From inspection of the diagrams 
that contain vertex corrections one sees that this entails at least three propagators with the 
mass of order T. Such diagrams are infrared finite. The diagram with a selfenergy insertion 
has the same infrared properties as the one without insertion so will be infrared finite h . In 
conclusion, no two loop kinetic term contributes 0(g 3 ) to the tension. 



6 The 0(g 3 ) contribution to the effective potential 

We have now eliminated the infrared divergencies from the potential up and including three 
loop order. By resumming the Z(N) invariant selfenergy H^(q) this procedure respects the 
Z(N) invariance of the potential. 

Although we have managed to get rid of the divergencies and keeping Z(N) invariance, 
there is a price to pay! Now the selfenergy is negative in a certain window of q values of order 
1, see fig. (6). So this window comes from integrating over hard momenta inside the selfenergy 
loop. The problem has not to do with the infrared scales. It is a generic problem occurring in 
quantum corrections to tunneling through a barrier 23 . However, for a window of loop strengths 
i < < | the mass stays positive over the whole range of q values. In particular in the large 
N limit we will be able to extract the potential without this problem. 

From the result for the masses of the Coulomb propagators in the previous section we find 
the contribution to the free energy. It will be a sum of terms of the type 

g 3 V 3 (q) = T I -i log(P + Ile(g)) (48) 
J ( 2n ) 

over all eigenvalues. 

We write for the ratio = r and introduce the function: 

/(r,g) = (l + 6rg(g-l))*-l (49) 

Note that this function is non-positive. The cubic potential, eq.(48), becomes in terms of this 
function: 

9 3 V 3 = -^-Tm 3 D N 2 (f(r, q)(l - rf + f(l - r, q)r 2 - ± (/( r , q ) + f(l - r ,q)- f(l, q))) (50) 

All terms proportional to cancelled out. The leading terms, proportional to ./V 2 , are well 
defined within the window | < r < | as seen from eq.(49). 



7 The tension to cubic order 

The cubic term in the tension is obtained from the minimization. The relevant expression was 
given in eq.(27): 

Pk = - [ dq(K(q)V(q))i (51) 
9 Jo 

h or will contribute <? 3 log(i) to V if q = 0(g 2 ) so 0(g 7 log(l/</) to the tension 
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n e (q) 



width of wall 




Loop momenta of order T 

Figure 6: The selfenergy for r = |. Note that the selfenergy is maximal for q = and q = 1 and equals there 
the (normalized) Debye mass. 



Combine the classical part of the kinetic term and the cubic part in V, eq.(50) and integrate 
over q. 

The tension becomes then to cubic order: 

Pr (T) = - 1.0341..(a s 7V) + (JL)!/(r)(a s AT)t + 0(a 2 s )) (52) 

The function 

is shown in fig. (7). It shows a small deviation from k-scaling, as it drops from the value 1.5 at 
r=0 to <~ 1.3 at r=0.5. As in the pressure 3 it contributes with a sign opposite to the 0(g 2 ) 
term. 

The attraction between loops becomes stronger due to the convexity of I(r). To see this, 
take the ratio of the tension of the k-loop and compare it to the k times the tension of the 
elementary loop with k=l: 

P±(T) 3 , , 1 

= (1 - r)(l + (^)HI(r) I(0))(a.N)* + 0(a 2 ) + 0{— 2 ) (54) 

The point is that the cubic correction has a negative coefficient, so that the ratio is smaller due 
to the presence of this correction. 
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Figure 7: The function I(r) in the result for the tension, eq.(52) 



8 Discussion 

We have computed the cubic correction to the tension of the 't Hooft loop. At first sight the 
occurrence of such a term seems counter-intuitive. One would expect the profile to provide an 
infrared cut-off on the order of \q\ = 0(T) well inside the wall. Only where the profile merges 
with the plasma ground state one expects the infrared divergencies to occur, as we discussed 
in section(4) for the next to leading correction to the Debye mass. 

This expectation is refuted by the presence of neutral gluonfields. They do propagate as 
massless fields in a constant background. But in the wall they can produce pairs of charged 
fields and interact therefore with the wall. The latter fact suggests that one should put the 
profile into the selfenergy, and then determine the free energy of the wall by computing the 
determinant of the operator. This method works not only for large N but for all SU(N) gauge 
theories. The mathematics of this method was used by Dashen, Hasslacher and Neveu 23 . ft 
can be worked out analytically and will be the subject of a future publication 24 . 

Concerning the magnetic components of the neutral gluons: do they get a self mass? As 
we show in the Appendix the magnetic selfenergy vanishes also in the presence of the wall. 
So a neutral magnetic gluon stays neutral, even in the presence of the wall. This means that 
their effects are only calculable by non-perturbative means. Or, what amounts to the same, to 
obtain a mass, one would need to put in a spatial Wilson loop. 

Still, it may be that all logarithmic effects of the type we encountered for the Debye mass 
can be captured by the mass induced into the off-diagonal magnetic glue. In other words the 
definition of the free energy in the bulk as a limit of our effective action where the Polyakov 
loop goes to its bulk value is not only manifestly gauge invariant but may also provide a well 
defined procedure to compute those logarithms. 

A last remark: our calculation applies to the domain walls discussed in ref. 25 . In five 
dimension the effects are g 4 log(l/g) instead of g 3 . 
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Appendix A: k-scaling of the two loop contribution 

In this appendix we prove the factor k(N — k) in front of eq.(35). The diagrams in fig. (2b) are 
all reduced by energy momentum conservation to the form: 

E /° lM /°' M E/iE/^ = E f aAt f aAt B2(C a )B 2 (C b ) (55) 

a,M J a J 6 a,M 

in an evident notation (d = 0, dj — d — Cj). Split the sum over the colour factor t into 
t = d, d diagonal and t = (i =/= j) . 

l)In the first sum a = ij, and b = ji, not diagonal. Then the completeness relation: 

N 1 

E f j ' ji ' d f d ' kl ' lk = -(S ik + 8n - 5 a - S jk ) (56) 

d=2 

says the first sum equals 

E r^r^Mc^MCji) = zhn- A)(s 2 (g)) 2 

ij,ji,d 

because on the path Y k (q) dj — ±q for 2k(N — k) combinations and because the sign of the 
argument dj is immaterial in the even function B^. 

Now the sum of t over ij. There are two different contributions. 

2a)First there is either a = d or b = d. Both give the same contribution 2k(N — fc)_B 2 (0)-B 2 ((z). 

2b)Then there is a = jl, and b = li. In case the indices i and j are in different sectors of Yfe 
either I equals one of the N — k indices i or one of the k indices j. In either case the summand 
equals B2(0)B2(q). And the coefficient in front is 2k(N — k)(N — 2), the factor N — 2 coming 
from the N — 2 ways we can put I. 

In case the the indices i and j are in the same sector we get 

(k(k - 1)(N -k) + (N-k)(N-k- l)k)B 2 (q) 2 
Adding up the contributions from 1), 2a and b) we get finally for the q dependent part: 

Y.\f a,b ' t \ 2 B2{d)B 2 (C b ) = Nk(N-k){B 2 {qf + 2B 2 (0)B 2 (q)) (57) 

a,b,t 

The contribution from the diagram in fig. (2c) equals: 

-4g 2 Nk(N -A0£i(?)£ 3 (<?) (58) 

where the first factor is the renormalization of the Polyakov loop and the second the derivative 
of the one loop potential. 

Added, they give the result in the text, eq.(34). 
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Appendix B: the selfenergy of the gluon in nontrivial background 

We want to prove eq.(42) and show that the magnetic part of the selfenergy vanishes at zero 
momentum in the presence of the wall. Look at fig. (5). The question is, how do wc identify the 
summand for fixed indices a and b in eq.(42)? For the first two graphs in fig. (5) this is rather 
obvious: the two internal lines carry the index a and b and will therefore have propagators 
w = 77, — ~ \ p> w ith C a = —Cb- For the tadpole graph we symmetrize the sole propagator 
in a and b: 

f c,a, bf c>,a, b J2 -1 = f Ca, bf c',a, b J2 I ( * + -L ) ( 59 ) 

In this formula the sum is over Matsubara frequencies only, not over colour indices. So the 
equality results from changing in the second term on the r.h.s. the summation from n into -n 
and using C a = —Cb- Then one finds: 



Use 



to find indeed: 



^C a ) = -T^j^lf-2ll)^ 2 (60) 
r d i p _ d r d i i 



T 2 
2 

The formula for the selfenergy in tensor form reads: 



(61) 



B 2 (C) = —B 2 {C) (62) 



n„„(C9 = -tJ2I Sfl« l a - 2 ^)i (63) 

n ■ ' a 

Contract the purely spatial part with 5ij and use eq.(61) to find it is identically zero, for any 
value of C and d. This means that the wall induces only a Z(N) invariant Debye mass on the 
Coulomb fields. 

Appendix C: eigenvalues of the selfenergy matrix 

We will diagonalize the selfenergy matrix, eq.(42) in the main text : 

n5(«) = ±g 2 f c - a - b B 2 (C a )f'' a - b (64) 

This matrix has no overlap between c = ij and d = d sectors as one can easily from the 
routing the colour indices in the self energy diagrams. So the matrix is divided into two sectors. 

The first sector is determined by the subset c = ij and d = ran defined by those phases 
Cij being identically zero on the tunneling path Y k (q). Recall that 

2-7T 

Cij = ir m{9))ii-{Y k {q)) jj ) 

With the definition Yk(q) = qYk we find therefore Cij = or ± 2irq. Here we are interested in 
ij such that dj = 0. There are 0(N 2 ) such pairs. 
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The second one is the one where the two propagators in the diagram are colour diagonal. 
It will beaiV— lbyTV— 1 matrix and will differ from one tunneling path to the other. The 
second sector involves only O(N) degrees of freedom and will be subdominant. 

Eigenvalues in the off-diagonal colour sector 

The one loop self energy n\j^ mn {q) is diagonal in ij and mn. So its value can be written, 
from eq.(42), as: 

U%\q) =2j29 2 \f ii ' jhM \ 2 T\B 2 (0) + B 2 (C hj )) (65) 
h 

Obviously after summing over h the Debye mass NB 2 (0) should appear at Chj = or 1. 
In the channel k one has two cases: 

i) The index pair (ij) is from the first N ~ k entries in the charge matrix Yfe. Then the running 
index h must be in the last k entries of , to give a nonzero contribution B 2 (Cih = B 2 (q). So 
£2(0) in eq.(65) is multiplied by N, the N possible values of the index h. But B 2 (dh) is non 
zero only for the k values of the index h. So the result is: 

U\f(q) = 2g 2 T 2 (NB 2 (0) + kB 2 {q)) = m 2 D (l + 6kB 2 (q)) (66) 
for any of the (N — k)(N — k — 1) pairs ij. 

ii) The index pair (ij) is from the last k entries in Y^. Then the index h must take on one of 
N — k values in the matrix Yk and the result becomes: 

H§\q) = 2g 2 T 2 (NB 2 (0) + (N - k)B 2 (q)) = m 2 D (l + 6(N - k)B 2 (qj) (67) 

for any of the k(k — 1) index pairs ij. 

Eigenvalues in the diagonal colour sector 



First a few general considerations. Of course any tunneling path differing from Yk(q) by a 
permutation P of the diagonal elements of Yfc(g) is physically the same: the potential should 
be identical. It is not hard to show that any such permutation amounts to an orthogonal 
transformation O acting on the orthogonal basis of the A^. So in particular the structure 
constants f d > % i>i % transform like: 

fd,P(ij),P(ji) = 0dA ,f'^^ (68) 

Now our selfenergy matrix reads: 

n K' (?) = ^9 2 f d ' ij ' ii B 2 (C ij )f d '' i " i (69) 

where we sum over the ij indices. 

It is then clear that the selfenergy transforms into OIl^O -1 , when computed from B 2 (Cp^ ) 
instead of B 2 (Cij). So its eigenvalues stay the same. And only the eigenvalues appear in the 
potential. Another property is that the eigenvalues are the same for charge conjugate paths 
Yk(q) and Yjv_fe(g), as expected. This comes about because a permutation changes one into 
the other up to a minus sign. This sign shows up in the argument of B 2 {dj) but is immaterial 
because B 2 is even. Like before it is useful to split B 2 into its value at q — and the rest: 

B 2 (q) = B 2 {{)) + B 2 (q) (70) 
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The first term contributes for every pair ij in the sum that constitutes the mass matrix. Hence 

2 

it comes in as proportional to the unit matrix with proportionality constant ^f- . The remaining 
term generates the non-trivial part of the selfenergy matrix: 

2 r' ,J '<rr'n-^rf' iJ ' (71) 

ij 

In contrast to the q-independent part this part of the mass matrix is in general non-diagonal. 
For k = 1 it is diagonal, for k = 2 it has off diagonal elements between d = N and d = N — 1, 
for k — 3 between d = N, N — 1 and N — 2. After diagonalization the non-trivial mass matrix 
has the eigenvalues (dropping the common factor Qm D B2{q): 

k (72) 

with multiplicity TV — k — 1. 
Then it has k — 1 eigenvalues 

(N - k) (73) 

and one eigenvalue 

N (74) 

As dictated by charge conjugation and periodicity modulo N the above set of eigenvalues 
does not change when we exchange k and N-k. 

A useful check on this result comes from taking the trace of the non-trivial selfenergy 
matrix, using ^2 d |J d »^>J J | 2 — 1. So this trace counts the number of nonzero eigenvalues in the 
adjoint representation of Yfc, which is 2k(N — k). So is the sum of the eigenvalues we found by 
explicit diagonalization of the selfenergy matrix above. 

Taking the results for the eigenvalues and multiplicities together one arrives at the set of 
eigenvalues and multiplicities in section(5). 
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